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Abstract 

A general theoretical framework is developed for inter- 
preting spectropolarimetric observations of optically thin 
emission line scattering from small dust particles. Spa- 
tially integrated and spatially resolved line profiles of both 
scattered intensity and polarization are calculated analyt- 
ically from a variety of simple kinematic models. These 
calculations will provide a foundation for further studies of 
emission line scattering from dust and electrons in such di- 
verse astrophysical environments as Herbig-Haro objects, 
symbiotic stars, starburst galaxies and active galactic nu- 
clei. 

Subject headings: ISM: dust, extinction — ISM: jets and 
outflows — ISM: reflection nebulae — polarization 

1 Introduction 

The Doppler shifts of optical emission lines which have 
been scattered by surrounding dust and electrons can pro- 
vide useful information about the kinematics, geometry, 
and physical conditions of astrophysical flows. In princi- 
ple, the scatterers can provide views of the line-emitting 
gas from different directions, allowing the 3-dimensional 
velocity of the emitting gas to be determined and reveal- 
ing sources which are obscured from direct view. Unfortu- 
nately, the interpretation of these scattered emission lines 
is not straightforward since, in general, the geometry of the 
scattering and the velocity of the scatterers is unknown. 
If spectropolarimetric observations are available, then de- 
termination of the relative orientation and velocity of the 
source and the scatterers becomes more plausible, since 
the scattered light will be partially polarized to a degree 
dependent on the angle of scattering and on the details of 
the scattering process. 

In this paper, a variety of very simple models of the 
scattering geometry and kinematics are investigated. Po- 
larimctric line profiles are calculated (both spatially in- 
tegrated and spatially resolved) for the cases of Rayleigh 
scattering by small dust particles, where the scattering 
phase function has a particularly simple form and it is fea- 
sible to derive analytic expressions for the scattered line 
profiles. 

Previous work on the problem of the scattering of emis- 
sion lines by dust or electrons has tended to involve de- 
tailed numerical radiative transfer modelling, often using 
Monte Carlo codes. Examples of the use of such tech- 
niques can be found applied to, e. g, electron scattering 
in Wolf-Rayet winds (Hillier 1991) or scattering in cir- 
cumstellar dust shells (Lefevre 1992; Meaburn, Walsh & 



Wolstencroft 1992). While such approaches have obvious 
merits, they are limited in their applicability, often need- 
ing to be specially tailored for each object, and it is diffi- 
cult to draw general conclusions from them. A simpler ap- 
proach is taken by Wood, Brown & Fox (1993), who model 
electron scattering in the circumstellar disks of Be stars, 
but, even here, numerical techniques are employed. In the 
present paper, numerical calculations are eschewed for the 
most part and scattered line profiles of intensity and po- 
larization are calculated analytically. The price that must 
be paid in order to make this approach tractable is to 
restrict one's attention to optically thin Rayleigh scatter- 
ing of narrow emission lines which originate from a point 
source. The chief advantage of this approach is that it 
makes it possible to investigate a wide range of models 
for the geometry and kinematics of the scatterers and to 
develop a broad understanding of the ways in which these 
affect the line profiles. This understanding can then pro- 
vide a firm foundation for further investigation into the 
effects of different phase functions (dust), thermal broad- 
ening (electrons), an extended source, absorption, multi- 
ple scattering and other complications. This development 
will be carried out in further papers of this series and 
detailed applications of the resultant techniques will be 
made to particular cases such as upstream dust scattering 
in Herbig-Haro objects (Noriega-Crespo, Calvet & Bohm 
1991), dust scattering in the jet of the symbiotic star R 
Aquarii (Solf 1992), dust scattering in the superwinds of 
starburst galaxies (Scarrott, Eaton & Axon 1991) and elec- 
tron scattering in active galactic nuclei (Miller, Goodrich 
& Matthews 1991). Some preliminary results in these ar- 
eas can be found in Henney (1992). 

2 Rayleigh Scattering 

Arbitrarily polarized light is fully described by means 
of the four-component Stokes vector (Stokes 1852) I = 
[/, Q, U, V] T , which can be functionally defined in terms 
of measurements of the light intensity as follows: / is the 
total intensity of the light; Q is the difference between 
the intensities after the light has passed through a linear 
dichroic plate with its transmitting axis aligned respec- 
tively perpendicular and parallel to a reference direction; 
U is the difference between the intensities with the same 
plate aligned at ±45° to the reference direction, and V is 
the difference between the intensities after the light has 
passed through a right and left handed circular dichroic 
plate respectively. Unpolarized light has Q = U = V = 
and fully polarized light has I 2 = Q 2 + U 2 + V 2 . Partially 
polarized light has a degree of polarization 

p=(Q 2 + U 2 + V 2 ) 1 / 2 /l (1) 
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and an angle of polarization x (with respect to the refer- 
ence direction in the definition of Q and U above) given 

by 



tan2x = C//Q. 



(2) 



Note that, since one must choose a reference direction cor- 
responding to x — 0, there is a certain arbitrariness in the 
specification of Q and U. In astronomy this reference di- 
rection is conventionally taken to be North (or sometimes, 
the galactic plane). The effect on the Stokes vector of 
a rotation of the reference direction by an angle <f> is to 
multiply it by the rotation matrix 



M(0): 
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cos 20 sin 20 

— sin 20 cos 20 

1 



(3) 



In the following discussion, it is imagined that light is 
emitted from a source at the origin and scattered by a sin- 
gle scatterer at [x, y, z] toward an observer who is in the 
far field and is looking down the negative z-axis (see Fig- 
ure [|). When the light is scattered, the components of the 
scattered Stokes vector are linearly related to those of the 
incident Stokes vector. Hence, the scattering cross-section 
can be written as a matrix. One way of representing the 
differential cross-section per solid angle ft is then 



3o-p 
16tT 



K, 



(4) 



where Co is the mean scattering cross-section and the ma- 
trix K (the scattering kernel; Pomraning 1973) contains 
the direction-dependent scattering information. The ma- 
trix <Tq is normalized so that 



<ToId£! = <7 I 



(5) 
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for any Stokes vector I. If Iq is the Stokes vector of the 
unscattered light at the source (in the direction of the scat- 
terer), then, by considering the solid angle at the source 
subtended by the scatterer, the scattered Stokes vector 
will be 



I s =<x o I /i? 2 



(6) 



where R 2 = x 2 +y 2 +z 2 . If the scatterers are much smaller 
than the wavelength of the light, then, so long as both are 
specified with respect to the scattering plane, the Stokes 
vector of the scattered light is related to that of the in- 
cident light via the Rayleigh scattering matrix (Chan- 



Figure 1: Geometry of single scattering. Light travels from a 
source at the origin (o) in a direction z' until it is scattered 
at a point b into a new direction z, towards the observer. 
Perpendicular axes x and y are in the plane of the sky such 
that yz is the observer's meridional plane. The perpendicular 
projection of B onto the z-axis is A and onto the iry-plane 
is C. The rectangle OABC thus lies in the scattering plane, 
containing the incident and scattered rays. The scattering 
angle between the two rays is O and the position angle of 
the scattering plane is ip. x' and y' complete an orthogonal 
set of axes with z' such that y'z' is the meridional plane for 
the incident light. The perpendicular projection of c onto 
the a/y'-plane is then D, which is coplanar with OABC, so 
that the position angle of the scattering plane in the primed 
coordinate system is ifr'. 
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drasekhar 1960, p 37): 



RO) = - 



2 e- 
2 e- 
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2 cos9 








2 cosG 



(7) 

where 6 is the scattering angle, the angle between the 
incident and the scattered ray. In general, if the scattering 
plane is at an angle ip to the reference direction, then R 
must be pre- and post-multiplied by appropriate rotation 
matrices M to give the scattering kernel K(O,"0). In our 
example, the scattering kernel is 



K(6, ip) =M(vr - ip)R(Q)M(-ip/) 



(8) 



where the angles ip and ip/ are shown in Figure |l| and 
the reference direction is taken to be the y-axis. If the 
incident light is unpolarized, then Q, U, and V are zero 
in the incident ray, so the M(— ip/) rotation need not be 
carried out and the Stokes vector of the scattered light is 
given by 

' ,T!l/o rX(e,^), 



AttR 2 ' 



where 



x(e,vo = : 



cos 2 e+i 

- cos 2tp sin 2 6 
sin 2ip sin 2 e 




(9) 



(10) 



The vector scattering phase function X(e,'0) is normal- 
ized so that, when integrated over all solid angles, the I- 
component becomes 47r and all other components become 
zero. Since single scattering of initially unpolarized light 
cannot produce any ^-component to the scattered Stokes 
vector, the ^-component will be dropped in the rest of the 
paper and a reduced Stokes vector I = [/, Q, U] T will be 
used. From equations (JTij) and (Q) it can be seen that the 
degree of polarization of Rayleigh scattered light is 



p = sin 2 e/(cos 2 e + l) 



(11) 



while the angle of polarization \ i s > using equation (g), 

X = 0±7r/2. (12) 

Hence the scattered light is always partially polarized per- 
pendicular to the plane of scattering and is fully polarized 
for scattering angles of 90° . 

3 Simple Scattering Geometries 



stationary dust cloud scattering light from a moving unpo- 
larized source at its center; second, that of an outflowing 
dusty wind; and third, that of a free-falling dusty inflow, 
the last two both scattering light from a stationary unpo- 
larized central source. In the first instance, the scattering 
cloud/wind/infall is assumed to be spherically symmetric 
and the light source is assumed to be a point, although 
the first of these assumptions is subsequently relaxed for 
the first two cases. 

Considering the spatially integrated light first, the three 
models are equivalent and it is evident on symmetry 
grounds that the scattered light will have no net polar- 
ization. This can be verified by taking the scattered in- 
tensity from a single scatterer (eq. ^) and integrating it 
over all the scatterers in the cloud. If the number density 
of scatterers is n(R) and the cloud radius is R c then this 
yields 



3ct /o 



n{R) 



16tt ./cloud R 2 



X(e,V)dF = r/ 



= rl 



(13) 

where r = a® J q c n dR is the scattering optical depth to the 
source. Of course, since t<1, the scattered light makes 
a negligible contribution to the total intensity unless the 
direct light from the source is obscured. 

3.1 Spatially Integrated Scattered Line 
Profiles 

3.1.1 Moving Source in a Stationary Scattering 
Cloud 

If it is now supposed that the source of light is moving 
with a velocity Uq with respect to the (stationary) dust 
cloud and that the light it emits is monochromatic (of 
frequency ^o), then the scattered light will be Doppler 
shifted by Ah* — (vq/c)R ■ Uq, where R is the unit vector 
in the direction of the scatterer from the source. One can 
define a dimensionless velocity shiftQ (frequency shift) of 
the scattered light as v = — R ■ Uq/uq = — cos 7. This is 
illustrated in Figure ^, in which the x-axis is taken to be 
along the projection of Uq onto the plane of the sky and 
the angle between Mo and the a;-axis is denoted by a. It 
is apparent that all scatterers that lie on the surface of a 
cone of opening half- angle 7 "see" the same Doppler shift 
(v = — cos 7) in the light from the source. This will be 
called the "isovelocity cone" and £ will signify the angle 



In this section, three simple classes of problems are con- 
sidered as examples of Rayleigh scattering: first, that of a 



1 In keeping with established convention, red shifts will be 
treated as positive. 
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around this cone, measured from a direction parallel to 
the y-axis (see Figure ||) . 

The spatially integrated scattered line shape can best 
be determined by considering the relationship between the 
two angular co-ordinate frames [7, £] and [O,^]. Equat- 
ing Cartesian components of a unit vector in the two co- 
ordinate frames gives 

x = sin sin ^ = sin 7 sin £ sin a + cos 7 cos a (14a) 
y = sin cos ip = sin 7 cos £ (14b) 
z = cos© = sin 7 sin £ cos a — cos7sino! (14c) 

so that, using equation (|l0|), the scattering phase function 
becomes 



X(7,0=4 



ao + ai sin£+a2 sin 2 £ 
b + bi sin£+fo 2 sin 2 £+63 cos 2 £ 
Ci sin £ + C4 sin £ cos £ 



(15) 



where 

a = sin 2 acos 2 7 + l; ai = — | sin 2a sin 27; a2 = cos 2 asin 2 7; 
bo = cos 2 a cos 2 7; 61 = | sin 2a sin 27; 



&2 =sin 2 a sin 2 7; 



63 = -sin 7; 



ci = — cos a sin 27; C4 = —2 sin a sin 7. 



To find the scattered intensity in a velocity range v± to V2 
it is necessary to integrate equation (||) over all the dust 
between the two isovelocity cones corresponding to v\ and 
V2- Hence, 



^.j p 7T — COS V2 t>ZTT 

I(vi->v 2 ) = -7^ / X(7,£)sin7d7d£. 

(16) 

Since di> = sin7d7, the specific intensity (per unit velocity 
range) is then 



Figure 2: Geometry of scattering from a moving source. The 
source velocity Uq lies in the xz-plane and is inclined at an 
angle a to the plane of the sky (xy). The radius from the 
source to the scatterer makes an angle 7 with Uq. The az- 
imuth of the scatterer, measured about Uq from a direction 
parallel to the y-axis, is £. The dotted lines show the "isove- 
locity cone", on which all scatterers see the same projected 
source velocity uq cos 7. 



I, = 



so that 



dl 


tI f 


dv 


^ Jo 



X(«,Od£ = -T/ 



a - 



7«2 



(17) 



I v = — T/ {2(l+t; 2 ) + (l-3i; 2 )cos 2 a}, (18a) 



— tIq (1 - 3v 2 ) cos 2 a, 



(18b) 



0. 



(18c) 



The polarization of the scattered light is then 
1 — 3v 2 cos 2 a 



Pv = 



2(l + u 2 ) + (l - 3v 2 )cos 2 a 



(19) 
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Figure 3: Line profiles of spatially integrated light, Rayleigh 
scattered by a stationary spherical dust cloud that surrounds 
a moving source, for various inclinations a of the source ve- 
locity to the plane of the sky. Solid line is specific intensity 
and dashed line is degree of polarization. 



while the angle of polarization \ 1S zero (perpendicular 
to the projected source velocity) when Q v is positive, or 
90° (parallel to the projected source velocity) when Q v 
is negative.^ These quantities are plotted in Figure || for 
various values of a. 

It should be noted that equation (|l^) is more general 
than equations (18a) and (18b) and applies to any distri- 



bution of scatterers that has cylindrical symmetry about 
the direction of the source velocity. This is because the 
number density n(R, 7, £) cancels out when taking the ra- 
tio of Q and I. (Cylindrical symmetry is still required, 
since it is assumed when performing the £ integral in 
eq. 0.) 

To illustrate how the intensity profile is modified if the 
scattering cloud is not spherical, the case of prolate and 
oblate ellipsoidal clouds will be considered. 



Ellipsoidal Scattering Clouds: If the scattering cloud 
is non-spherical, then the chance of a source photon being 
scattered is dependent on its initial direction. For a pro- 
late ellipsoid of eccentricity e, with symmetry axis aligned 
with the source velocity, this can be translated into a di- 



2 Notice that the polarization is here allowed to take the sign 
of Q, rather than being constrained to be positive. This is a 
useful convention in instances like this one, where \ can only 
take one of two values, but can be distracting otherwise. 



Figure 4: Same as Figure g but for prolate spheroidal scat- 
tering clouds. Different line types correspond to different 
eccentricities, as indicated in the key. 



(20) 



rection dependent optical depth of 

t( 7 ) _/ l-e 2 sin 2 cA 1/2 
\ 1 — e 2 cos 2 7 / ' 

where r* is the optical depth along the line of sight from 
the observer (corresponding to "/ — tt/2 — a). Hence, the 
intensity of the scattered light is given by 

1/2 



1 2-2 

1 — e sin a 
1 — e 2 v 2 
,:2\ 



x {2(l + u 2 ) + (l-3v 2 )cos 2 a}. (21) 

A similar treatment for an oblate ellipsoid yields 

1/2 



1 — e 2 sin 2 a 



(1 - e 2 ) + e 2 v 2 , 
x {2 (l + u 2 ) + (l -3u 2 ) cos 2 a} . (22) 

The resultant line shapes are plotted for various values of 
e and a in Figures ^ and |5|. The polarization is not shown 
since it is the same as for a spherical cloud. It is apparent 
that prolate ellipsoids produce line shapes in which most 
of the intensity is concentrated toward v = ±1 whereas 
oblate ellipsoids produce line shapes strongly peaked at 
v = 0. The reason for this can be seen if one considers 
the distributions of scatterers in the two cases: in pro- 
late clouds, most of the scatterers are concentrated along 
the direction of the source velocity, either "upstream" or 
"downstream", while in oblate clouds they are concen- 
trated toward the plane perpendicular to the source ve- 
locity and hence see little Doppler shift. It will be noticed 
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Figure 5: Same as Figure fi| but for oblate spheroidal scatter- 
ing clouds. 



that the integrated intensity (with respect to v) is not the 
same in each case. This is because the intensities are nor- 
malized with respect to t*Iq, where r* is the scattering 
optical depth to the source along the line of sight. Gen- 
erally, t* is not equal to the angle-averaged optical depth 
(weighted by the phase function) f , although it is this lat- 
ter that determines the integrated scattered intensity. For 
example, a prolate ellipsoid at a = or an oblate ellipsoid 
at a~7r/2 would each have r* < f. 

3.1.2 Constant Velocity Scattering Wind 

If the scattering wind is supposed to flow radially outward 
at constant speed u w , then the Doppler frequency shift in 
the scattered light will be Ai/ — (vo / 'c)u w (l+R ■ z), where 
z is the unit vector along the z-axis (from the source to- 
ward the observer). In the same manner as for the moving 
source model, a dimensionless velocity shift v = 1+R ■ z is 
introduced which again defines isovelocity cones. In this 
model, however, the opening half-angle of the cone is the 
scattering angle and this is related to the velocity shift 
by cos O = v— 1. Determination of the scattered line profile 
is hence far simpler than in the preceding section, giving 



- 2v- 





(23) 



The scattered line is unpolarized because the isovelocity 
cones project onto the plane of the sky as circles. This line 
shape is illustrated in Figure || and it can be seen that it 
is identical to the moving source model with a = 7r/2, 



Figure 6: Spatially integrated line profile from a constant- 
velocity spherically symmetric wind of scatterers. 



except for a shift in the velocity origin. If the wind is not 
spherically symmetric but is conical or disk-like in form, 
then the symmetry of the isovelocity cones is broken and 
polarized line profiles can result. 

Outflowing Bicone: The outflow axis of the wind is 
taken to make an angle a with the plane of the sky and 
the wind density is assumed to be independent of direction 
within the two cones of opening half-angle S c about this 
axis, and zero outside the cones. These scattering cones 
will intersect a given isovelocity cone along 0, 2 or 4 radii 
with position angles (^i . . . Vl^) given by 



sin VP i 



sin ^2 



sin^q = sin\I/4 = 



cos 8 C — cos 9 sin a 
sin 9 cos a 

- cos 5 C — cos 9 sin a 
sin 9 cos a 



(24) 



where 9 = cos 1 \v — 1|. The scattered Stokes vector is 
therefore given by 




#2 /•*4 

X(u,*)d*+ / X(w,*)d* 



which can be integrated to give the following: 
Iv = -^-tI (v 2 ~2v + 1)F 1 

l07T 



Qv - —tI {2v-v 2 )F 2 

107T 



u v = 0, 



(25) 

(26a) 
(26b) 
(26c) 
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where 



F, = < 



2tt 

2tt-2(*i- 




[if a-6>n/2-S] 

■* 3 ) [if 0-a > n/2-5] 

( [if a+9 > tt/2 + 6 ( 27a ) 
\ or a+(9 < tt/2-(5] 

[otherwise] 



Fa 





sin2* x - 


sin 2*i 



sin 2* 



3 



[if a-6>TT/2-5] 

[if 0-a > n/2-5] 

f [if a + 9 > ir/2 + 5 ( 27b ) 

\ or a + 9 < k/2-5] 

[otherwise] 



In Figures and& examples of these line shapes are given 
for narrow (<5 C = 0.2 radians) and wide (<5 C = 1.0 radians) 
cones. 

Outflowing Disk: If the wind is confined to a thin equa- 
torial disk, of angular half-thickness Sd then the scattered 
line profiles can be calculated by subtracting the result 
for a biconical wind from that of a spherical one so that 
5 C \ = tt /2 — 5 c . Neglecting terms higher than quadratic in 
<$d (thin disk) then gives 



T 3 . . v 2 -2v + 2 



{2v-v 2 



■sin a 



1/2 



(28a) 



3 (l+sin 2 a)(2u-i; 2 )-2sin 2 a 
8 cos z a(2u — ir)(2i> — ir — sin a) 1 '^ 

where a is now the angle between the normal to the disk 
and the plane of the sky. Examples of these line shapes 
are given in Figure ||. 

3.1.3 Free-falling Inflow 

If a dusty flow is assumed to fall radially inwards from 
rest at infinity towards a gravitating body of mass M* 
then, neglecting any deceleration, the infall velocity will 
be given by u in = (i? /-R) 1/2 uo wher e u = {2GM*/ Rq) 1 / 2 
is the velocity at the inner cut off radius Rq (which may be 
identified with the grain destruction radius). The inflow 
is supposed to have an outer radius R c and mass con- 
servation dictates that the dust number density have the 
form n = ria(Ra/R) 3 / 2 . The Doppler frequency shift of the 
scattered light is Al> = — (uo/c)ui n (l+R ■ z), in a similar 
manner to the outflowing wind but with the opposite sign. 



Figure 7: Integrated spectra of intensity (solid line) and de- 
gree of polarization (dashed line) for a scattering biconical 
wind with cone opening half-angle S c = 0.2 and at various 
inclinations a of the cone axis to the plane of the sky. 
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Figure 8: Integrated spectra of intensity (solid line) and de- 
gree of polarization (dashed line) for a scattering biconical 
wind with cone opening half-angle 5 C = 1.0 and at various 
inclinations a of the cone axis to the plane of the sky. 



Figure 9: Integrated spectra of intensity (solid line) and de- 
gree of polarization (dashed line) for a scattering wind in the 
form of a thin disk at various angles a between the normal 
to the disk and the plane of the sky. Note that the intensity 
diverges at v = 1 ± cos a and so the values at these points 
are a function of the velocity resolution used in the plots. 
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Figure 11: Spatially integrated scattered line profiles from a 
spherically symmetric free-falling inflow with ratios of inner 
to outer flow radius eg = 0.05, 0.25 and 0.5. 



Figure 10: The isovelocity scattering surfaces for a free- 
falling inflow of dust are surfaces of revolution formed by 
rotating the above curves about the line of sight to the cen- 
ter of the flow. The curves are marked with the dimensionless 
Doppler shift induced in emission lines scattered by dust lying 
on them. In this example, the ratio of inner to outer dust 
radius is eo = 0.2. As can be seen, the isovelocity surfaces 

1/2 

are closed for v < — 2e n ~— 0.89. 



from the point of view of the observer. Example line pro- 
files are shown in Figure |ll] for various values of eo- It 
can be seen that for high eo (corresponding to a small 
spread in dust velocities with radius) the profile is similar 
to that from the constant velocity wind but with the sign 
of the Doppler shift reversed. As eo is decreased there is 
a greater proportion of slow-moving dust and so the line 
profile becomes increasingly skewed toward v = 0. 



In this case however the dust velocity is not constant, so 
the dimensionless velocity shift is dependent on the dust 
radius: v = — (e /^) 1 / 2 (l + cos6), where R = R/R C and 
e a = R /R c - This means that the isovelocity surfaces are 
no longer cones but have the more complex shapes shown 
in Figure [l0|. By integrating the scattered intensity over 
these isovelocity surfaces, it is possible to determine the 
spatially integrated scattered line profile of the infall as 



if v < -2e, 



1/2 



8(l- £o /2 ) 



G 2 if v > -2e 



1/2 



where 



G 2 



= lne -2,(e 1 / 2 -l)-^ 2 ( 



.2 



2 In ( --v 



1 9 

-2 + 2v + -v 2 



(29) 

(30a) 
(30b) 



Note that the scattered line profiles are unpolarized be- 
cause the isovelocity surfaces are circularly symmetric 



3.2 Spatially Resolved Scattered Line 
Profiles: Position- Velocity Diagrams 

In this section, an analytic approach is used to calculate 
the Stokes intensities of the Rayleigh scattered light re- 
solved both spatially and in velocity for two of the classes 
of models presented above. 

3.2.1 Moving Source in a Stationary Scattering 
Cloud 

For the purposes of this section, the scattering cloud will 
be taken to be spherical and homogeneous. The observer 
is looking down the negative z-axis (as in Figure ^) and 
the position of a scatterer in the cloud is characterized by 
dimensionless coordinates [x,y, z] where x — x/R c etc. . 
Straightforward geometry then shows that the dimension- 
less Doppler shift of light scattered at [x, y, z] is 



x cos a +z sin a 



(x 2 +y 2 + z 2 ) 1/2 ' 



(31) 
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Hence, solving this equation for z, a line of sight [x,y] 
intersects a given isovelocity cone v in zero, one or two 
places given by 

{z a , Zb} = I x sin a cos a ± v [x 2 (l — u 2 ) 



+y 2 (sin 2 a — u 2 )] ^ 2 ^ x (y 2 — sin 2 a) 



so long as v satisfies 



(32) 



icosai [l-(.T 2 +y 2 )] 1/2 < \v\ < ( 1- 

.i — ii i 

(33) 

where the positive sign in the equation and inequalities 
applies to z a and the negative sign to §b- The left-hand 
side of equation (|3^) reflects the fact that most lines of 
sight will not intersect the source velocity vector, while 
the right-hand side is the result of the finite size of the 
scattering cloud. The scattered line profile along a line of 
sight can then be calculated as 



y 2 cos 2 a 



-y 



1/2 



I, = 



dl 

dv 



Att9 2 



E 



X(a:,y, Zj) 
x 2 +y 2 + z 2 



dzi 



dv 



(34) 



where Fo is the flux of the source and 9 C is the apparent 
angular radius of the scattering cloud. The phase function 
in Cartesian coordinates can be written as 



X(x,y, z)-. 



x 2 +y 2 + 2z 2 

t2 ^.2 



4(x 2 + y 2 + z 2 ) 
and from equation ( |3l"l ) it follows that 



x- - y 
-2yz 



(35) 



dz 




dv 





(x 2 + y 2 + z 2 ) 



3/2 



(x 2 +y 2 ) sin a — zxcosa 



(36) 



It is then necessary to eliminate z from equation ( |34| ) us- 
ing equation (^). For the general case this leads to very 
complicated expressions which are best calculated by com- 
puter but results for two special cases are presented here. 

Source Velocity in Plane of the Sky (a — 0): If 

the dimcnsionlcss velocity shift v lies between x and 



xf (x 2 +y 2 ) 



1/2 



then 



Iv = 



3tF [x 2 (2 - v 2 ) - v 2 f] 
8n8 2 x[x 2 {l-v 2 )-v 2 y 2 ] 1/2: 



3tF v 2 (x 2 -y 2 ) 



iir9 2 x [x 2 (1 — v 2 ) — v 2 y 



(37a) 
(37b) 



U v = 



-3tFq v 2 xy 



Ane 2 x[x 2 (1 - v 2 ) - v 2 y 2 ] 1/2 ' 



(37c) 



otherwise /„ = Q v = U v = 0. The degree and angle of 
polarization are given by 



v 2 (x 2 + y 2 ) 
\2x 2 - v 2 (x 2 + y 2 )\' 



Pv 



X = 2+tan- H - 



(38a) 



(38b) 



Note that x is independent of v, which is true for all values 
of a. 

Source Velocity Along Line of Sight (a = 7r/2): For 
a = tt/2 the scattering geometry has circular symmetry 
from the point of view of the observer so it is convenient 
to introduce a dimcnsionlcss impact parameter r — (x 2 + 
y 2 ) 1 ' 2 . Then, if |v| < (l-r 2 ) 1 /2 j 



F, = 



3tF (1+v 2 ) 

3TF {x 2 -y 2 )(l-v 2 )^ 2 
167r6> 2 r 3 ' 

-6rF xy(l -v 2 ) 1 ' 2 



(39a) 



(39b) 



(39c) 



167r6» 2 r 3 

otherwise = Q v = U V = 0. The degree of polarization is 

1-v 2 



Pv = 



1 + V 2 



(40) 



and the angle of polarization is still given by equa- 
tion (38b). Note that the polarization is independent of 
position on the cloud for a given velocity, as is the in- 
tensity (apart from a scale factor). However, the range of 
allowed velocities is strongly dependent on the position on 
the cloud. 

Example Spectra and Position- Velocity Diagrams: 
In Figure [l2| the positions on the scattering cloud of six 
apertures used in constructing example spectra, three slits 
and three point apertures, are shown. In Figures 13 
to [l8], the spectra and position-velocity diagrams from 
these apertures are shown for different values of a. These 
spectra include the broadening effects of the frequency 
profile of the source emission line and of atmospheric see- 
ing. Both are assumed to be Gaussian in form; the source 
profile is taken to have a FWHM of 0.16ito and the seeing 
profile to have a FWHM of O.110 c . 
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Figure 12: The positions, on a spherical cloud, of the syn- 
thetic apertures used in § 3.2.1. Two-dimensional position- 



velocity diagrams are presented for Apertures A-C (Fig- 
ures |l3|-|l5|) and one-dimensional spectra for Apertures D-F 
(Figures p^6[-|l8|). 



3.2.2 Constant Velocity Scattering Wind 



Spherically symmetric wind: For a constant velocity 
wind (assuming a constant dust-gas ratio), mass conser- 
vation requires the dust number density to be of the form 



(41) 



where no is the number density at an inner cut-off radius 
Ro = eoR c - The scattering optical depth to the source is 
then given by 

T = <T n R c eo(l-e ) (42) 



so that the the line profiles are given by 
t tF e X(x,y,z) 



4tt9 2 l-e (r 2 + z 2 ) 2 



dz 
dv 



(43) 



Note that in this instance a line of sight can only intersect 
an isovelocity cone in, at most, one place and that the 



situation is very similar to the a = it/2 case of § 3.2.1 
Hence, if (l-(r/e ) 2 ) 1/2 < < (1-r 2 ) 1 / 2 then 



h = 



Pv 



3TF e (2v-v 2 ) 1 / 2 (v 2 -2v + 2) 



167r6» 2 (l-e )r 3 



2v-v 2 
v 2 -2v + 2' 



(44a) 
(44b) 



Figure 13: Doppler- and seeing-broadened position-velocity 
diagrams for a stationary cloud, scattering an emission line 
originating in a moving source. Aperture A: centrally placed 
slit, aligned with the source velocity. Dimensionless Doppler 
width w = 0.1. Dimensionless seeing width a s = 0.066. Inten- 
sity contours are logarithmic, with successive contours corre- 
sponding to a ratio of 2 1 / 2 . Value of lowest contour is given 
on each diagram in units of tIq and is 2~ 6 times the peak 
value. Also shown are vectors whose lengths are proportional 
to the degree of polarization and whose orientation shows the 
position angle of the polarization. 



12 



Figure 14: Same as Figure 

parallel to the source velocity. Lowest contour is 2 
the peak value. 



13 but for Aperture B: offset slit, Figure 15: Same as Figure [l3] but for Aperture C: offset slit, 
4 5 times perpendicular to the source velocity. Lowest contour is 2~ 4 - 5 
times the peak value. 
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Figure 16: Doppler- and seeing-broadened spectra for a sta- 
tionary cloud, scattering an emission line originating in a 
moving source. Aperture D: centered on the source. Di- 
mensionless Doppler width w = 0.1. Dimensionless seeing 
width a s = 0.066. 



Figure 18: Same as Figure [l6| but for Aperture F: centered 
to the side of the source. 



otherwise I v = (vis still given by eq. [38b] ) . As in the 
a = ir/2 case of § 3.2.1, variation in the line profiles with 
position is chiefly caused by the conditions for intersection 
between the line of sight and isovelocity cone. 



Outflowing Bicone: With a conical wind, the intensity 
and polarization profiles will still be given by equation (Q) 
but the range of allowed velocities for the scattered light 
from a given line of sight is subject to the additional con- 
dition that the line of sight must intersect the cone of 
scatterers. To simplify matters, only lines of sight that 
intersect the cone symmetry axis (y = 0) are considered. 
In this case, the conditions that there be scattered flux at 
velocity v from a line of sight x are given in Table 1, where 
all three conditions must be satisfied using either entirely 
the left terms in curly brackets or entirely the right. The 
position velocity diagrams so obtained (corresponding to a 
slit placed along the projected axis of the conical outflow) 
are illustrated in Figures |l^ and ^ for both a narrow 
cone and a wide cone at various inclinations to the plane 
of the sky. 



Outflowing Disk: Since the disk is assumed thin, the 
spatially resolved scattered line profile at each position on 
the disk will just be a single spike. Hence, the position- 
Figure 17: Same as Figure |l6| but for Aperture E: centered velocity diagram for a slit placed along the projected minor 
upstream of the source. axis of the disk (y = 0) will be a ridge at v = l — cosa for 

negative x and a ridge at u = l+cosa for positive x, both 
with a degree of polarization sin 2 a/(l + cos 2 a), whereas 
for a slit along the projected major axis of the disk (x = 0) 
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TABLE 1 

Conditions for there to be Scattered Flux from an Outflowing Bicone at a Point x, v in 

Position- Velocity Space 



(l-(5/6 ) 2 ) 1/2 < < (1-x 2 ) 

max [cos(a+<5), 0] < {x,—x} < cos (min[a — S, 0]) 

sin(a-*) [*{*,-*}><>] | f s H a + 5) ^j^^ 

> < {1— y, < < 

sin(a + <5) [if{z,-2}<0] J [ 1 [otherwise] 



Note. — See text for explanation of the symbols. 



the position-velocity diagram is merely a ridge at v = 
with 100% polarization. 

4 Discussion 

The calculations presented in this paper have been se- 
lected because they are amenable to analytic treatment 
and hence represent extreme idealizations of situations 
likely to be encountered in astrophysical objects. Nev- 
ertheless, they demonstrate in a simple fashion the type 
of scattered line shapes that will be produced in different 
situations. 



means that a source moving toward the observer would be 
seen to have an unpolarized red wing of width 2uq and 
a source moving away from the observer would have a 
similar unpolarized blue wing, while a source moving in 
the plane of the sky would have polarized red and blue 
wings. In this latter case, the polarization is highest at 
the extremes of the wings, drops to zero and then rises 
again toward line center (although in the center dilution 
by the intrinsic light will reduce the observed polarization 
substantially) . Departure from spherical symmetry of the 
scattering cloud (Figures | and §) causes the blue- and 
red-shifted scattered light to become more (less) intense 
than the unshifted light if the cloud is prolate (oblate). 



4.1 Moving Source 



The moving source model (§§ |3.1.1| and 3.2.1) is applicable 
to any case in which a line-emitting plasma is moving with 
respect to a dusty environment, such as Herbig-Haro ob- 
jects. This is discussed in much greater detail in Henney, 
Raga & Axon (1994) and Henney, Axon & Raga (1994), 
Papers II and III of this series. It may also furnish an 
alternative explanation to that proffered in Solf (1992) for 
the enormous line widths observed in knots of the colli- 
matcd outflow from R Aquarii. 

It is found that the spatially integrated line shapes from 
these models show scattered wings extending a distance 
equal to the source velocity uq to both sides of the rest 
frequency of the line (Figure ||) . These wings are polarized 
to a degree dependent on the inclination of the source 
velocity to the line of sight, being highest when the source 
is moving in the plane of the sky and zero if the source 
is moving directly toward or away from the observer. Of 
course, the Doppler shift of the direct light from the source 



Of the results for spatially resolved line shapes that are 
presented in § 3.2.1 , it is perhaps those for Aperture A of 
Figure [l2] (corresponding to a narrow slit placed along a 
diameter of the cloud, parallel to the direction of motion 
of the source) that are the most interesting (Figure |l^) . 
For the source velocity in the plane of the sky, these 
show a characteristic "double triangle" morphology to the 
position- velocity diagram of the scattered light. The scat- 
tered light is blue-shifted to the right (upstream) of the 
source and red-shifted to the left (downstream). The po- 
larization is highest for the scattered light that undergoes 
the largest Doppler shift (red or blue) since it is the dust 
directly in front of or behind the source that scatters light 
through 90°, leading to maximum polarization (Eq. 11 1). 
As the angle a between the source direction and the plane 
of the sky increases, some upstream dust appears to the 
left of the source and some downstream dust to the right. 
Eventually, when the source is moving directly toward or 
away from the observer, this leads to the situation illus- 
trated in the bottom-right panel of Figure O in which the 
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Figure 19: Scattered position velocity diagrams from a slit 
placed along the axis of a narrow outflowing bicone of open- 
ing half-angle <5 C = 0.2 radians at various inclinations a to the 
plane of the sky. 



Figure 20: Same as Figure 19 but for a wide outflowing 



bicone of opening half-angle <5 C = 1 . radians. 
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position- velocity diagram is symmetrical about x — 0. In 
this latter case, it is the scattered light which undergoes 
no Doppler shift which has the highest polarization since 
the dust in the plane of the sky is now directly to the sides 
of the motion of the source. Note that, for both a = and 
a = 7r/2, the area of the position-velocity diagram in which 
the direct light from the source will lie (x = 0, v = for 
a = 0, u = ±l for a = ±7r/2) is well away from the area of 
highest polarization so there will be little dilution of the 
scattered light there. Note also that, unlike in the spa- 
tially integrated case, the scattered light is polarized for 
all inclinations of the source velocity. 

4.2 Scattering Wind 



These models (§§ 3.1.2 and [3.2. 2| ) are applicable to any 



case in which chromospheric emission lines are scattered 
from a stellar wind. This may be dust scattering, as is 
possible in evolved late-type stars such as Mira variables 
(Anandarao, Pottasch & Vaidya 1993; Lefevre 1992), or 
electron scattering, e. g. in Be stars (Wood, Brown & Fox 
1993), B[e] supergiants (Boyd & Marlborough 1991), or 
Wolf-Rayet stars (Schulte-Ladbeck et al. 1992). They are 
also relevant to the scattering by dust and electrons in out- 
flows from active galactic nuclei (Miller & Goodrich 1990). 
Of course, for electron scattering the thermal broadening 
of the electron velocity distribution will significantly mod- 
ify the results presented here for all but the highest Mach 
number winds (the thermal speed of the electrons will ex- 
ceed the bulk wind velocity so long as the Mach number 
M < 2(toh/to c ) 1 ^ 2 ~80), a fact that is ignored in Wood, 
Brown & Fox (1993). The extension of these results to 
Thomson scattering from thermal electrons will be pre- 
sented in a forthcoming paper. 

One obvious result is that some asymmetry of the wind 
is necessary for the integrated scattered line profile to be 
polarized. For a spherical wind (Figure ||), the scatter- 
ing merely produces an unpolarized red wing, extending 
up to twice the wind velocity u w . If the wind is concen- 
trated toward the polar directions (Figures [?] and |8|), then, 
for outflows in the plane of the sky, one finds a polarized 
red bump at a redshift of u w whose width depends on the 
degree of collimation of the wind. As the outflow axis 
is rotated toward the observer (increasing a), the bump 
splits into two, which move apart to redshifts of zero and 
2u w as a — > it/2 and whose polarizations diminish to zero. 
For cquatorially enhanced winds (approximated as a disk; 
Figure ^|), the intensity profiles of the scattered lines for 
a given inclination of the symmetry axis a can be seen to 
be qualitatively similar to the profiles for a narrow cone 
(Figure f7j) with inclination ir/2—a, as would be expected. 
The polarization profiles do not follow this pattern, how- 



ever, and show a maximum at a red shift of u w and a 
general decrease in polarization as the disk changes from 
an edge-on to a face-on orientation. 

The spatially resolved position-velocity diagrams are 
not illustrated for the spherically symmetric wind but they 
are very similar to those for the moving source model with 
a=n/2 (bottom-right panels of Figures H H), except for 
three differences: 

1. The Doppler shifts will range from to 2u w instead 
of from — u w to u w . 

2. In the spectrogram from Aperture A, there will be 
an elliptical "hole" in the center because of the inner 
cut-off to the density distribution. 

3. The scattered brightness will fall off more rapidly with 
distance from the center of the cloud because of the 
inverse-square density distribution in the wind. 

With the conical wind (Figures [l9| and |2(]) , the position- 
velocity diagrams are the same as for the spherical wind 
but with certain portions masked out. The sizes and po- 
sitions of the non-empty regions of position- velocity space 
depend on the opening angle and orientation of the cones. 



4.3 Scattering Inflow 



The free-falling inflow model (§ 3.1.3) is applicable to 
the scattering haloes of young stellar objects, which of- 
ten show density profiles indicative of an infalling envelope 
(Weintraub et al. 1992). 

The spatially integrated line shapes (Figure |ll|) can be 
seen to depend sensitively on the value assumed for the 
ratio of inner to outer dust radius eo (see discussion after 
cq. ) . It should be remembered that the dimensionless 
Doppler shift v is scaled to the infall velocity at the inner 

— 1/2 

cut-off radius, which is proportional to e , so that the 
velocity scales of Figure [ll] are different if all the inflows 
are assumed to have the same outer radius and central 
mass. 

Calculations of the spatially resolved line profiles from 
this model are not presented here because they are too in- 
volved for the simple analytic approach used in this paper. 
In addition, it is likely that many young stellar objects 
have optically thick haloes (Whitney & Hartmann 1993) 
and in such cases the single scattering model presented 
here is not strictly applicable. These issues will be ad- 
dressed in Paper III, in which Monte Carlo simulations of 
multiple scattering will be made and where the effects of a 
non-zero angular momentum for the infall will be treated. 



17 



5 Summary 

This paper has presented a small collection of line pro- 
files (both intensity and polarization) that result from the 
optically thin Rayleigh scattering of an emission line by 
dust in its environment. The calculations have all been 
performed analytically and this has necessarily restricted 
the complexity of the models employed. Nonetheless, they 
form a basis for understanding the ways in which the ge- 
ometric and kinematic relations between the line source 
and the scatterers determine the scattered line profiles. 

Although, the calculations have been performed for the 
optically thin case, for one of the kinematic models con- 
sidered (moving source in a stationary cloud) large op- 
tical depths will make very little difference to the scat- 
tered line profiles. This is because there is no relative 
motion between the scatterers themselves and hence the 
only Doppler shift is that induced by the first scattering. 
Of course, for the wind and infall models, the optically 
thick case will be quite different, with much more extended 
wings to the profiles due to multiple scattering. However, 
the polarization will also be reduced compared with the 
single scattering case so that the optically thin models are 
perhaps more likely to be relevant to sources with observ- 
able polarization changes across their line profiles. 

One effect that has been ignored in this paper and that 
can have a significant effect on the line profiles is that 
of a non-Rayleigh scattering phase function for the dust. 
This is probably important at optical wavelengths, where 
the scattering phase function can be quite forward-peaked 
unless the dust grains are smaller than usual. This issue 
will considered in depth in Paper III. 
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